Void Ellipticity Distribution as a Probe of Cosmology 
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Cosmic voids refer to the large empty regions in the universe with a very low number density 
of galaxies. Voids are likely to be severely disturbed by the tidal effect from the surrounding dark 
matter. We derive a completely analytic model for the void ellipticity distribution from physical 
principles. We use the spatial distribution of galaxies in a void as a measure of its shape, tracking 
the trajectory of the void galaxies under the influence of the tidal field using the Lagrangian per- 
turbation theory. Our model implies that the void ellipticity distribution depends sensitively on the 
cosmological parameters. Testing our model against the high-resolution Millennium Run simulation, 
we find excellent quantitative agreements of the analytic predictions with the numerical results. 
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PACS numbers: 98.65.Dx, 95.75.-z, 98.80.Es 



Introduction. — One of the most fundamental goals in 
physical cosmology is to understand the large-scale struc- 
ture in the universe. The organization of the large-scale 
structure begins at the scale of galaxies and appears to 
follow a hierarchical pattern up to the scale of superclus- 
ters 

The study of the large-scale structure has been so far 
focused mostly on the bound objects like clusters and 
superclusters. However, more recent observations indi- 
cate that the universe is in fact a collection of bubble- 
like voids separated by sheets and filaments, at the dense 
nodes of which the clusters and superclusters are located 
Q. Hence, to address the large scale structure in the 
universe, it seems essential to account for the presence 
and characterization of voids. 

According to the standard gravitational instability the- 
ory, voids originate from the local minima of the initial 
density field and expand faster than the rest of the uni- 
verse. The extremely low-density of the observed voids 
[d v ~ —0.9) supports this scenario Q. A common expec- 
tation based on this standard scenario is that voids are 
likely to have quite spherical shapes 0, 0, H, Hi- How- 
ever, a recent systematic analysis of simulation data has 
revealed that the shapes of voids are in fact far from 
spherical symmetry 0, [1] . 

Confronted with this somewhat unexpected result, 
Ref. Q claimed that the gravitational tidal field be re- 
sponsible for the nonspherical shapes of voids. They ex- 
plained that since the voids have very low-density, they 
must be more easily disturbed by the tidal effect from 
the surrounding matter. 

Reference [9J has investigated the tidal effect on voids 
analytically using the linear tidal torque theory, which 
was inspired by [8[. It is shown by them that the tidal 
field indeed has a significant influence on voids, gener- 
ating non-radial motions of matter that make up voids. 
They tested their analytic model against high-resolution 
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N-body simulation and found remarkably good agree- 
ments between the analytical and the numerical results. 

The success of the analytic model of [§] has two crucial 
implications. First, the non-radial motions of galaxies in 
voids generated by the tidal effect would cause the non- 
sphericity in void shapes, as hinted first by ||. Second, 
since the Lagrangian theory has turn out to work well in 
predicting the tidal effect on voids, it may be also possible 
to model the nonspherical shapes using the Lagrangian 
theory that depends only on the initial conditions. 

One difficulty in modeling the nonspherical shapes of 
voids, however, lies in the fact that there is no well de- 
fined way to determine the shape of a void since a void is 
an unbound system having no clear-cut boundary. Nev- 
ertheless, given that our purpose is not to construct the 
full geometrical shape of a void but to quantify the degree 
of the deviation of its shape from spherical symmetry, a 
practical way that we choose here is to use the spatial 
distribution of the void galaxies as a measure of the non- 
spherical shape of a void. If a void has a spherical shape, 
then the spatial distribution of its galaxies will show a 
more or less isotropic pattern. If it has a nonspherical 
shape, the void galaxy spatial distribution will also de- 
viate from the isotropic pattern. In this practical way, 
the nonsphericity of the shape of a given void may be 
quantified by the ellipticity of the spatial distribution of 
the galaxies that it contains. 

In this Letter we aim for constructing a complete an- 
alytic model for the void ellipticity using the Zel'dovich 
approximation (the first order Lagrangian perturbation 
theory) and investigate if the void ellipticity distribution 
can put constraints on the cosmological parameters. 

The Analytic Model — Assuming that the density field 
is smoothed on the void scale, we use the Zel'dovich ap- 
proximation 

.10] to describe the trajectory of a galaxy lo- 
cated in a void region as x = q— D(t)V v I'(q). Here x and 
q represent the Eulerian and the Lagrangian positions of 
a void galaxy, respectively, D(t) is the linear growth fac- 
tor and \l/(q) is the perturbation potential smoothed on 
the Lagrangian void scale, Rl- Note that here a galaxy 
ac.kRi treated as a particle moving under the influence of 'J. 
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Applying the mass conservation law, pc? 3 x = pd 3 q to 
the Zel'dovich approximation, the mass density of a void 
region, p v , can be written as 



Pv 



p 



(1-A 1 )(1-A 2 )(1-A 3 ) : 



(1) 



where Ai,A2,A3 (with Ai > A2 > A3) are the three 
eigenvalues of the tidal tensor, Ty, defined as the sec- 
ond derivative of The dimensionless residual overden- 
sity of the void region, S v , equals the sum of the three 
eigenvalues: S v = Y^=i A*. 

Eq. ([!} implies that under the influence of the tidal 
field, Tij, the void region will experience a triaxial ex- 
pansion. Let pi,p2,pz (with pi > P2 > P3) represent the 
three semi-axis lengths of the inertia momentum tensor 
of the galaxies that make up void , and let [i and v be the 
two axial ratios defined as p = pilpx and v = ps/pi. The 
ellipticity of a given void region in terms of the minor- 
to-major axis ratio as e = 1 — v. 

According to Eq. (TT]), the two axial ratios are related 
to the three eigenvalues of (Tij) as 
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The distribution of e can be found from the distribution 
of u, which in turn can be derived from the conditional 
distribution of the three eigenvalues of the tidal tensor 
under the constraint of 6 V = Y)i—i Ai- 

In Ref . (llj . the unconditional joint probability density 
distribution of {Ai, A2, A3} was derived. On the void scale 
Rl, it is written as 



p(Ai, A 2 ,A 3 ;cr fl ,J 
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2< 



where K x = Ai + A 2 + A 3 , K 2 = AiA 2 + A 2 A 3 + A3A1, 
and <7ji l is the rms density fluctuation related to the 
dimensionless linear power spectrum A 2 (k) as 
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A 2 {k)W 2 (kR L )d\nk, 
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where W(kRL) represents a smoothing function with the 
Lagrangian filtering radius of Rl- Here, we adopt the 
approximation formula given by [l2| for A 2 (k) and the 
top-hat spherical filter for W(kRL). 

By Eq. ([3]), the conditional joint probability density 
distribution of the void axial ratios is found as 



p(fi,u\S v ;(7 RL ) = A p[\ 1 (p,is),\ 2 (p,is)\6 v ;a RL 
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Here p(X±, X 2 \5 V ; ct Rl ) is the conditional probability den- 
sity of Ai and A 2 provided that a given void has an av- 
erage density S v on the scale of Rl- This conditional 



probability has the following form: 
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where the functional forms of X\(fi,v) and A 2 (/z, v) are 
given as [HI 
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and the normalization constant, A, satisfies the con- 
straint of A J p(n, v\8 = 5 v ;a RL )dpdv = 1 

The distribution, p(v) is now straightforwardly evalu- 
ated by integrating Eq. ([5]) over /i as 



p{p,v\5 = 6 V ; o-R L )dfj,. 



(9) 



Then, the distribution of the void ellipticity, p(s), is noth- 
ing but p(e) = p(l — v). Fig. [1] plots p(e) at two different 
Lagrangian scales: Rl = 5 and 8/i -1 Mpc in the upper 
and the lower panel, respectively. Here we assume a A- 
dominated cold dark matter cosmology (£1a = 1 — Cl m ) 
and consider four different sets of the key cosmologi- 
cal parameters: Q m (matter density parameter), h (di- 
mensionless Hubble parameter), as (amplitude), and n s 
(power spectrum slope). 

As can be seen, the ellipticity distribution depends 
quite sensitively on the cosmological parameters. Note 
that for the cases of low f2 m and high h, voids tend to 
have more spherical shapes, which in fact can be under- 
stood as follows. The void ellipticity distribution is an 
outcome of the counterbalance between the tidal effect of 
the dark matter and the expansion of the universe: The 
tidal effect disturbs the void shapes from spherical sym- 
metry, while the expansion of the universe resists such 
disturbance.lt is also worth noting that the void elliptic- 
ity distribution moves toward the low ellipticity as the 
smoothing scale increases. It indicates that the larger a 
void is, the more spherical it is. 

Numerical Test. — To examine the validity of our an- 
alytic model, we test its prediction against the Millen- 
nium Run simulation for which the cosmological param- 
eters are specified as fi m = 0.25, n s = 1, <r 8 = 0.9 and 
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In our previous work Q we have already identified 
voids in the Millennium Run catalog using the Hoyle- 
Vogeley (HV02) void-finder algorithm [ll|. For the de- 
tailed description of the void-finding scheme and the 
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FIG. 1: Probability density distributions of the void ellipticity 
predicted by our analytic model for the four different sets of 
cosmological parameters at two different smoothing scales: 
Rl = 5 and 8ft~ 1 Mpc in the upper and the lower panel, 
respectively. For each case, the void density contrast is set at 
8 V = -0.9. 



properties of the identified voids, see [9j. Among these 
voids, we select only those voids which contain more than 
5 galaxies. A total of 23652 voids are found to contain 
more than 5 galaxies. 

For each selected void, we first calculate the pseudo 
inertia momentum tensor, (iy), defined as 

1 



hj — 



Nva 



^ ; X ai X aj' 
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where N vg is the number of galaxies belonging to a given 
void, and {x 9 ai ) is the position vector of the a-th void 
galaxy. As mentioned in Q, the inertia momentum ten- 
sor is computed without weighting the position vectors 
with mass since we are interested in the geometry of 
voids. 

Diagonalizing {hj), we find the three eigenvalues, ii, 
h, and /3 (with h > 1% > h), which are related to 
the semi axis- lengths of a given void as h cx (j>\ +p%), 
h oc {pi+pD, I3 oc {j>2+p\) Hi- Now, the void ellipticity, 
e = 1 — pa/pi, can be written in terms of h,h, h as 
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Measuring e of each selected void using Eq. (fTT|) , we de- 
termine the void ellipticity distribution. Since the distri- 
bution depends on the size of void, we divide the whole 



TABLE I: The range of the effective radius (i?jg m and 
Kg ), the mean effective radius (-RpJ, the mean density 
(<5v), the mean ellipticity (e), and the mean number of the 
galaxies {Nvg) of the voids belonging to the four bins. 
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FIG. 2: Illustration of the four examples of the voids identified 
in the Millennium Run simulation. Each panel corresponds 
to the four bins. The shaded regions and the dots represent 
the areas of the voids and the positions of the void galaxies 
in them projected onto the two dimensional plane. 



sample of voids into the four bins according to their effec- 
tive radius, Re, which is defined as Re = (3ry/47r) 1 / 3 
where Ty is the total volume of a void measured by the 
Monte-Carlo method devised by [lij]. Each bin is chosen 
to include equal number of voids. The range of Re, the 
mean values of Re, 8 v , and s for the four bins are listed 
in Table Q] As shown, e decreases as the mean R e in- 
creases. That is, the larger a void is, the more spherical 
shape it tends to have. 

Fig. [2] illustrates the projected images of the four sam- 
ple voids selected from the four bins, respectively. In 
each panel the dots represent the positions of the void 
galaxies and the shaded regions correspond to the areas 
of the voids. Fig. [3] plots the void ellipticity distribu- 
tions for the four bins. In each panel, the histogram and 
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FIG. 3: The ellipticity distributions of the voids belonging 
to the four bins. In each panel the histogram represents the 
numerical data points from the Millennium Run simulation 
while the solid line corresponds to the analytic prediction. For 
the analytic results, the values of the Lagrangian scale radius, 
Rl, are evaluated by the relation of Rl = (1 + Sv) 1 ^ 3 Re- 



the solid line represent the numerical and the analyti- 
cal results, respectively. For the analytic distribution, 
the Lagrangian scale radius, Rl, is given analytically as 
Rl = (1 + SvY^Re- As can be seen the analytic and 
the numerical results are in excellent quantitative agree- 
ments for all cases, proving the validity of our model. It 
is worth emphasizing here that our model is completely 
analytical without having any fitting parameter derived 
from principles. It depends only on the initial cosmolog- 
ical parameters. 

Discussion. — The ellipticity distribution of galaxy 
clusters has been used recently to constrain cosmology 
1^ . 17 . ljj ]. For the case of clusters, however, the elliptic- 
ity distribution cannot be evaluated analytically from the 
initial conditions because the nonlinear clustering pro- 
cess tends to modify the cluster ellipticities significantly 
in the subsequent evolution. Although [13[ attempted to 
derive an analytic fitting model for the cluster ellipticity 
distribution, their model was shown to be not in good 
quantitative agreements with the numerical results. 

In contrast, we have for the first time shown here that 
the void ellipticity distribution can be evaluated fully an- 
alytically using physical principles and thus be a better 
indicator of the initial conditions of the universe since the 
void ellipticities are less vulnerable to nonlinear modifi- 
cations. By applying our analytic model to real data that 



will be available from future large galaxy surveys , it may 
be possible to constrain the cosmological parameters pre- 
cisely in an independently way. 

It is, however, worth nothing that the Hoyle-Vogeley 
algorithm works very well if a void is found as a spherical 
low-density region. Nevertheless, wc expect that our nu- 
merical result would not change sensitively with respect 
to the void-finding algorithm, given the fact that the el- 
lipticity of a void here is defined in terms of the spatial 
distribution of galaxies that make up the void but not by 
the boundary shape of a void. 

The Millennium Run simulation used in this paper was 
carried out by the Virgo Supercomputing Consortium 
at the Computing Centre of the Max-Planck Society in 
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comment. This work was supported by the New Faculty 
Settlement fund of the Seoul National University. 
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